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We show that the large band offsets between GaN and InN and the heavy carrier effec- 
tive masses preclude the use of the Virtual Crystal Approximation to describe the electronic 
structure of Gai-^In^N/GaN heterostructures while this approximation works very well for the 
Gai-^In^As/ GaAs heterostructures . 
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The use of semiconductor alloys proves to be neces- 
sary to adjust the electronic parameters to the design 
of specific devices. Gai_. r In c N alloys inserted between 
GaN would be ideally suited to cover the entire near 
infrared/ultra- violet spectrum Despite the lack for 
translation invariance, one may very often describe the 
electronic behavior of an alloy by means of a band struc- 
ture, the so-called Virtual Crystal Approximation Q 
(VCA). Close to it, the Coherent Potential Approxima- 
tion allows to introduce part of the disorder effects in the 
VCA scheme, which leads to a damping of the VCA Bloch 
states. An example of alloys which is well described by 
the VCA/CPA is Gai-Jn^As. The VCA has also been 
used to describe the electronic properties of Ga(In)N het- 
erostructures 

In this letter, we shall show that the Gai_ x In x N sys- 
tem behaves in a radically different way. A numerical 
computation of the electronic states of Gai-^In^N/GaN 
quantum wells and quantum dots will show very large 
differences from the VCA predictions. Instead, we shall 
show that Gai-^In^As/GaAs electronic states are very 
close to the VCA predictions. The difference between 
the two alloy families arises from the much larger band 
offsets between GaN and InN than those between GaAs 
and InAs and from the heavier effective masses in the ni- 
tride system than in the arsenide system. We note that 
the large potential offset between In and Ga has led Kent 
and Zunger 5] to predict that nitride based alloys cannot 
be described by models which neglect fluctuations. 

In the following, we use one band effective mass Hamil- 
tonians to describe the electron and hole kinematics. 
The effective mass is slightly anisotropic in the conduc- 
tion band (m ez — 0.184 mo, m e / / = 0.166 mo) and 
anisotropic in the valence band (m^ z = 1.1 mo, m^n — 
0.5045 m ) HQ- The conduction AE C (valence AE V ) 
band offsets between GaN and InN is taken equal to 
1.8 eV (0.9 eV) while those between GaAs and InAs are 
taken equal to 0.41 eV and 0.29 eV respectively ||. Thus, 
there is a much larger energy fluctuation in the nitrides 
than in the arsenides when in a given unit cell a Ga atom 
replaces an In atom. The two heterostructures we shall 
be dealing with are a 3.2 nm thick Gao.83Ino.17N/GaN 
quantum well and a Gao.83Ino.17N truncated pyramid 
with hexagonal basis (6 nm side), 2.6 nm height and 



a basis angle of 30° embedded into a GaN matrix and 
floating on a 1.1 nm thick wetting layer (see Fig.QJ). A 
pyramid with the same geometrical parameters made of 
Gao.5Ino.5As and embedded into GaAs will be consid- 
ered for comparison (50% of In instead of 17% was con- 
sidered to get bound states to the pyramid). It is well 
known that nitride heterostructures contain huge inter- 
nal fields 0, 0- This will be modeled by assuming 
that there exist piecewise constant electric fields which 
are oriented along the Z direction and have a magnitude 
of 2.45 MV/cm in Ga sslno 17N and -0.1 MV/cm in GaN 
6,7]. 




Figure 1: Left: VCA band diagram and wave functions in 
quantum well in the growth direction. Right: geometry of a 
quantum dot with its wetting layer 

The quantum well problem will be reduced to a two 
dimensional (2D) calculation by considering an effective 
in - plane potential which is the average of the actual 3D 
potential weig hted by \$>(Z)\ where * (Z) is the enve- 
lope function solution of the ID problem in the presence 
of the electric field and a square well confining poten- 
tial with a depth of 0.31 eV and 0.15 eV (= 0.17A.E C 
or 0.17A^„) (see Fig.JTJ). This effective 2D potential is 
then supplemented by a hard box potential which con- 
fines the carrier in the layer plane into a square with 
dimensions L x L = 36 nm x 36 nm. A given sample is 
then randomly generated by filling each unit cell either 
by In or by Ga corresponding to a potential energy of ei- 
ther (In) or AE C (Ga) for electrons and either (In) or 
AE V (Ga) for holes. For the quantum dots a complete 3D 
diagonalization was undertaken enclosing the pyramid in 
a box with a square basis with a side of 14.9 nm and a 
height of 7.9 nm. The box is then subdivided into smaller 
boxes of atomic size where the potential is supposed to be 
constant. The wave functions are then developed on the 



2 



sinus solutions leading to a Hilbert space of dimension 
75000. 

For the QDs, two situations have been considered: ei- 
ther uncorrelated sites where the probability that a given 
cell is occupied by In is x ( = 17%) irrespective of the oc- 
cupancy of the other sites or site correlations which favor 
the In cluster formation by allowing a probability 

p > x for a cell to be occupied by an In if its first neigh- 
borings are already occupied by at least one In. The only 
samples retained in the analysis are those such that the 
In fraction in the whole QD lies in the interval [0.17, 0.18]. 

The first few eigenstates of those 2D (40 states) and 
3D (6 states) problems are evaluated by means of exact 
diagonalizations using the Lanczos algorithm. Once the 
eigenstates are known one can calculate various averages 
like the mean position, the mean square deviation to it, 
etc.. 

While the properties of a single sample have virtually 
no chance to be compared to any experiment, the aver- 
age properties, obtained by taking the arithmetic mean 
of the outcomes of a given physical quantity in a series of 
samples, have more physical substance since they can (in 
principle) be compared to the results obtained on ensem- 
bles of quantum dots or quantum wells. This is for these 
average properties that one can discuss the effectiveness 
or failure of the VCA. One important quantity is Z? e , the 
average electron density of states which is : 
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where N is the number of samples (which differ from one 
another by the locations of the In atoms) while e e i(j) is 
the i th eigenvalue of the j th sample and T has been taken 
equal to 1 meV. Here we replace the usual delta peaks by 
Lorentzian distributions for clarity. Equally important 
(and measurable) is the average electron - hole optical 
density of states defined as : 
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In the finite QW structure the VCA results in peaks re- 
lated to the confinement box which fulfill: 
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where Vb = xAE c is the VCA average potential, e\ z is the 
lowest confining energy in the Z direction and L =36 ran. 

Fig.J2J shows for the QW structure the mean square 
deviation 5a plotted versus the mean location (X) for the 
electron (upper panel) and the hole (lower panel) along 
the X axis (the results are equivalent for the Y axis) for 
50 samples, together with the VCA results. The In dis- 
tribution is at random (x = 0.17). It is obvious that (X) 
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Figure 2: Calculated mean square deviations 8 a plotted ver- 
sus mean carrier position (X) for the VCA (a), electrons (b) 
and holes (c) in a finite Gao.83Ino.17N/GaN quantum well. 
Random alloys and 50 samples. 
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Figure 3: Optical density of states averaged over N = 50 
Gao.83Ino.17N/GaN quantum wells. Random alloys. Dashed 
lines: VCA peaks. 



and 5a for holes differ strongly from the VCA predic- 
tions while the electrons with their lighter mass are less 
readily confined by In rich fluctuations and thus display 
5a ((X)) which are closer from the VCA results than the 
holes. Fig.||3J shows the electron - hole optical density of 
states D e h(e) versus the energy e averaged over N — 50 
samples. D e h is proportional to the light absorption co- 
efficient if we neglect electron-hole interaction, which is 
justified here due to the strong e-h separation by the elec- 
tric field (770 meV compared to the Coulomb interaction 
energy of the order of 44 meV). The VCA results in a 
rounded (by T) staircase - like D e h for an infinite QW 
while the large box produces VCA peaks whose position 
follows eq. Js3 ■ We see that the random alloy results in a 
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broadening of these peaks which is quite important. Note 
for instance the significant (D e h > 0.05) band tail which 
develops down to « 15 meV below the nominal edge of 
the ground peak. This bandtail corresponds to increas- 
ingly localized states around In clusters, in particular for 
the holes. The high energy decrease of the optical density 
of states is unphysical. It reflects the high energy cut-off 
of the electron and hole eigenvalues. Notice that the pho- 
toluminescence line is found experimentally at ~2.6 eV 
in Gao.s3Ino.17N/GaN with which our results are co- 
herent. 
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Figure 4: Electron density of states averaged over N = 100 
Gao.83Ino.17N/GaN pyramids, (n) is the mean number of first 
In neighbors for an In. From top to bottom: VCA results 
(divided by 6), random alloy (n) = 1.0, and increasingly seg- 
regated alloys (n) = 1.2; (n) = 1.6; (n) = 2.1. The energy 
zero is taken at the barrier monolayer next to the bottom of 
the wetting layer. 
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Fig.QJ shows the calculated electron density of states 
for Gao.83Ino.17N/GaN pyramid (100 samples); the opti- 
cal density of states displays similar results although the 
hole states can be unbound in the pyramid. The VCA 
results shows a twofold degenerate excited state located 
some 90 meV above the ground state. It lies some 80 meV 
below a narrow triplet. The random alloy case shows 
both a redshift and a considerable broadening. Note in 
particular the tendency towards the closing of the gaps 
between the remnants of the VCA peaks. This effect is 
more pronounced when In segregation is introduced in 
the calculations. In addition, increasing the segregation 
produces densities of states which bear less and less re- 
semblance with the VCA result. In contrast (fig.©) the 
Gao.5Ino.5As/GaAs pyramids (50 samples) display den- 
sities of states which are very close to the VCA results. 
Even though the 50% alloy should display the maximum 
of disorder, we see that the calculated DOS features are 
very well separated from each others, and much narrower 
than the ones found in the 17% nitride alloy. Note in 
particular that the bandtail extends only a few meV's 
below the VCA peaks. Thus, if one knows the shape of 
a Gao.5Ino.5As pyramid, fig. I© shows that it does make 
sense to attempt to identify the peaks, to fit their energy 
difference etc. On the other hand, the huge deviations 
from the VCA results in the 17% nitride dot demonstrate 
that any peak attribution in this system is rather elusive. 
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Figure 5: Electron density of states averaged over N = 50 
Gao.5Ino.5As/GaAs pyramids. Solid lines: random alloy. 
Dashed lines:VCA results. 



